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ABSTRACT 



Recently, Dixon et al. ( [1998 ) have re-analyzed the EGRET data, finding a statistically 



significant diffuse 7-ray emission from the galactic halo. We show that this emission can 
naturally be explained within a previously-proposed model for baryonic dark matter, in which 
7-rays are produced through the interaction of high-energy cosmic-ray protons with cold H2 
clouds clumped into dark clusters - these dark clusters supposedly populate the outer galactic 
halo and can show up in microlensing observations. Our estimate for the halo 7-ray flux turns 
out to be in remarkably good agreement with the discovery by Dixon et al. ( 199§| ). We also 
address future prospects to test our predictions. 
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1. Introduction and outlook 



As is well known, the galactic halo chiefly consists of dark matter. A natural possibility - repeatedly 



considered in the past (Silk 1993, Carr 1994) - is that baryonic dark matter makes a substantial contribution. 



A few years ago, we recognized that the Fall & Rees theory for the formation of globular clusters (Fall 
& Rees 1985, Kang et al. 1990, Vietri & Pesce 1995| ) also leads to the existence of dark clusters - made 
of brown dwarfs [j] and cold self-gravitating clouds - at galactocentric distances R > 10 kpc (De Paolis et 
1995a| - |l9954 |l998a| ) 0. Accordingly, the er halo is populated by globular clusters whereas the outer 



al 

halo is dominated by dark clusters. Contrary to the case of globular clusters, a large amount of residual 
gas should remain clumped into the dark clusters, as brown dwarfs fail to generate the strong stellar winds 
which expel the leftover gas from globular clusters. Moreover, although the clouds under consideration are 
mainly made of H2, we expect them to be surrounded by an atomic layer and a photo-ionized "skin" (De 
Paolis et al. 1998a). We stress that the presence of cold self-gravitating clouds in the halo is a characteristic 



feature of the model in question. Remarkably enough, quite recently it has been pointed out (Walker 



& Wardle 1998) that cold self-gravitating clouds of the considered kind naturally explain the "extreme 



scattering events" associated with compact radio quasars (Fiedler et al. g_987j) . 

Our proposal was also motivated by the discovery of microlensing events towards the LMC (Alcock et 



al. 1993, 1997, Aubourg et al. 1992). The first-year data were manifestly consistent with the assumption 



that MACHOs are brown dwarfs even within the standard (isothermal) galactic model. Unfortunately, the 
present situation is much less clear. An option is that the halo resembles more closely a maximal disk 
rather than an isothermal sphere, in which case MACHOs can still be brown dwarfs ^ (see also Binncy 
1998| ). A more intriguing possibility has recently been suggested by Kerins & Evans (1998) and naturally 
fits within our model. As the initial mass function evidently changes with the galactic distance R, it can 
well happen that brown dwarfs dominate the halo mass density without however dominating the optical 
depth for microlensing 



A few months ago, Dixon et al. (1998) have re-analyzed the EGRET data concerning the diffuse 7-ray 
flux with a wavelet-based technique. After subtraction of the isotropic extragalactic component and of the 
expected contribution from the Milky Way, they find a statistically significant diffuse emission from the 
galactic halo. At high-galactic latitude, the integrated halo flux above 1 GeV turns out to be ~ 10~ 7 — 10~ 6 



7 cm 2 s 1 sr 1 , which is slightly less than the diffuse extragalactic flux (Sreekumar et al. 1998 ) 



Our aim is to show that the diffuse 7-ray emission from the galactic halo discovered by Dixon et al. 



(1998) can naturally be explained within the considered model. Basically, the idea is that cosmic-ray (CR) 
protons in the galactic halo scatter on halo clouds, thereby producing the observed 7-ray flux. 



1 Although we concentrate our attention on brown dwarfs, it should be mentioned that red dwarfs as well can be 
accommodated within the considered scenario. 



2 Similar mod els ha ve also been proposed by Ashman & Carr (1988), Ashman (1990), Fabian and Nulsen (1994, 1997), 
Gerhard & Silk (1996) and Kerins ( 19974 Ujg7t ). 

3 It should be kept in mind that a large fracti on of M ACHOs (up to ~ 50% in mass) can be binary systems, thereby counting 
as twice more massive objects (De Paolis et al. 1998b). 

4 Notice that also the considered clouds can contribute to microlensing events (Draine 1998). 
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2. Cosmic ray confinement in the galactic halo 

Unfortunately, neither theory nor observation allow nowadays to make sharp statements about the 
propagation of CRs in the galactic halo 0. Therefore, the only possibility to get some insight into this issue 
rests upon the extrapolation from the knowledge of CR propagation in the disk. Actually, this strategy 
looks sensible, since the leading effect is CR scattering on inhomogeneities of the magnetic field over 



scales from 10 2 pc down to less than 1CP 6 pc (Berezinsky et al. 1990 ) and - according to our model - 
inhomogeneities of this kind are expected to be present in the halo, because of the existence of gas clouds - 
with a photo-ionized "skin" - clumped into dark clusters []. 

As is well known, CRs up to energies of ~ 10 6 GeV are confined in the galactic disk for ~ 10 7 yr. CRs 
escaping from the disk will further diffuse in the galactic halo, where they can be retained for a long time, 
owing to the scattering on the above-mentioned small inhomogeneities of the halo magnetic field PL 

Indirect evidence that CRs are in fact trapped in a low-density halo has recently been reported. For 



example, Simpson and Connell (199S) argue that, based on measurements of isotopic abundances of the 
cosmic ratio 26 A1/ 27 A1, the CR lifetimes are perhaps a factor of four larger than previously thought, thereby 
implying that CRs traverse an average density smaller than that of the galactic disk. 



A straightforward extension of the diffusion model (Berezinsky et al. 1990| ) implies that the CR escape 



time T es ^f from the halo (of size Rh much larger than the disk half-thickness) is given by 



h % (1) 

csc ~ 3D H (E) ' 1 ' 



where Dh{E) is the diffusion coefficient. 



We recall that - for CR propagation in the disk - the diffusion coefficient is D(E) ~ Do (E/7 GeV) ' 3 
cm 2 s _1 in the ultra-relativistic regime, whereas it reads D(E) ~ Dq ~ 3 x 10 28 cm 2 s _1 in the 



non-relativistic regime (see Berezinsky et al. 1990 ). As a matter of fact, radio observations in clusters of 
galaxies yield for the corresponding diffusion constant D a value similar to that found in the galactic disk 
(Schlickeiser, Sievers & Thiemann 1987) [| So, it looks plausible that a similar value for Dq also holds on 



intermediate scale lengths, namely within the galactic halo. In the lack of any further information on the 
energy-dependence of Dh(E), we assume the same dependence as that established for the disk. Hence, 
from eq. ([!]) - with Rh ~ 100 kpc - we find that for energies E < 10 3 GeV the escape time of CRs from the 
halo is greater than the age of the Galaxy i ~ 10 10 yr (notice that below the ultra-relativistic regime T e ^ 
gets even longer). As a consequence - since the CR flux scales like E~ 2 - 7 (see next Section) - protons with 
E < 10 3 GeV turn out to give the leading contribution to the CR flux. 



5 We stress that - contrary to the practice used in the CR community - by halo we mean the (almost) spherical galactic 
component which extends beyond ~ 10 kpc. 

6 Indeed, typica l values of the dark cluster radius are ~ 10 pc, whereas typical values of the cloud radius are ~ 10" 5 pc (De 
Paolis et al. 



1998e) 



7 A similar idea has been proposed with a somewhat different motivation by Wdowczyk & Wolfendale (199!:) 



Moreover, we note that averag e mag netic field values in galactic halos are expected to be close to those of galaxy clusters, 



i.e. in the range 0.1 - 1 fiG (Hillas 1984) 
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We are now in position to evaluate the CR energy density in the halo, getting 



Pen 



4:7rR 3 H 



0.12 eV cm" 



(2) 



where Lq ~ 10 41 erg s 1 is the galactic CR luminosity. Notice, for comparison, that p^ R turns out to be 
about one tenth of the disk value (Gaisser 1990] ). 

We remark that we have taken specific realistic values for the various parameters entering the above 
equations, in order to make a quantitative estimate. However, somewhat different values can be used. For 
instance, Rr may range up to ~ 200 kpc (Bahcall, Lubin & Dorman 1995), whereas Dq might be as large 



as 



10 2 



consistently with our assumptions. Moreover, Lq can be as large as 3 x 10 erg s 



(Volk, Aharonian & Breitschwerdt 1996 ). It is easy to see that these variations do not substantially affect 
our previous conclusions. 



3. Gamma-ray emission from halo clouds 

We proceed to estimate the total 7-ray flux produced by halo clouds clumped into the dark clusters 
through the interaction with high-energy CR protons. CR protons scatter on cloud protons giving rise 
(in particular) to pions, which subsequently decay into photons. We expect negligible high-energy (> 100 
MeV) 7-ray photon absorption outside the clouds, since the mean free path is orders of magnitudes larger 
than the halo size. 

An essential ingredient is the knowledge of both p^ and the CR spectrum Q^ R (E) in the galactic 
halo. According to the discussion in the previous Section, we take Pq R ~ 0.12 eV cm -3 . As far as &q r (E) 
is concerned, we adopt the following power-law 

^cr(E) - ^ (Sv) P articles cm " 2 s_1 sr_1 » ( 3 ) 

where the constant A is fixed by the requirement that the integrated energy flux (in the range 
1 GeV < E < 10 3 GeV) agrees with the above value of Pq R - The choice of a is nontrivial. As an orientation, 
the observed spectrum of primary CRs on Earth entails a ~ 2.7. However, this conclusion cannot be 
extrapolated to an arbitrary region in the halo (and in the disk), since a crucially depends on the diffusion 
processes undergone by CRs. For instance, the best fit to EGRET data in the disk towards the galactic 
centre yields a ~ 2.45 (Mori 1997] ), thereby showing that a increases as a consequence of diffusion. In the 



lack of any direct information, we conservatively take a ~ 2.7 even in the halo, but in the Table we report 
for comparison some results for different values of a. As can be seen, the flux does not vary substantially. 

Let us next turn our attention to the evaluation of the 7-ray flux produced in halo clouds through the 
reactions pp — > 7r° — > 77. The source function q 7 (> E^,p,l,b) - yielding the photon number density at 
distance p from Earth in the direction (Z, b) with energy > E 1 - is 

47r f 00 

q 7 (> E 7 ,p,l,b) = pH 2 (p,l,b) dE p <&(? R (E p ) a in {pi ab ) <n 7 {E p ) > 7 cm -3 s _1 , (4) 

m p Je p (E-,) 

where the lower integration limit E p {E 1 ) is the minimal proton energy necessary to produce a photon with 
energy > E 7 , <Ji n {pi a b) is the inelastic pion production cross-section, n 7 (E p ) is the photon multiplicity ^] 



For the inclusive cross-section of the reaction pp — > tt° — > 77 we employ the parameterization given by Dermer (1986) 
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and Ph 2 (p, Kb) is the halo gas density profile which reads [[]] 

PH 2 {x,y,z) = f Po (q) - T - — =- — " , (5) 

a z + x z + y z + (z/q) z 

for \f x 2 + y 2 + z 2 /q 2 > R m i n , otherwise it vanishes. Here R m in — 10 kpc is the minimal galactocentric 
distance of the dark clusters in the galactic halo, / denotes the fraction of halo dark matter in the form 
of gas, po(q) is the local dark matter density, a = 5.6 kpc is the core radius and q parametrizes the halo 
flattening. For the standard spherical halo model po(q — 1) ~ 0.3 GeV cm~ 3 , whereas it turns out that e.g. 
p (q = 0.5) ~ 0.6 GeV cur 3 . 

Because dV = p 2 dpdfl, it follows that the 7-ray flux per unit solid angle produced in halo clouds and 
observed on Earth from the direction (Z, b) is 

1 rP2(l,b) 

$ ° M (> E~,l,b) = — / dpqJ>E 1 ,p,l,b) 7 cm" 2 s" 1 sr" 1 , (6) 

47r -Wb) 

where typical values of p\ and pi are 10 kpc and 100 kpc, respectively. 



4. Discussion and conclusions 

Our main result - which follows directly from eq. (^|) - are maps for the intensity distribution of the 
7-ray emission from baryonic dark matter in the halo. In order to make the discussion definite, we take 
0.5. 

In Figure 1 we show the contour plots in the first quadrant of the sky (0° < Z < 180°, 0° < b < 90°) for 
the 7-ray flux at energy E 1 > 1 GeV $ 7 DM (> 1 GeV). Corresponding contour plots for E^ > 0.1 GeV are 
identical, up to an overall constant factor equal to 8.74 (again, this follows from eq. (Q)). 

Figure la refers to a spherical halo, whereas Figure lb pertains to a q = 0.5 flattened halo. We see 
that - regardless of the adopted value for q - at high-galactic latitude ( E' 7 DM (> 1 GeV) lies in the range 
~ 6 — 8 x 10~ 7 7 cm -2 s _1 sr -1 . However, the shape of the contour lines strongly depends on the flatness 
parameter. Indeed, for q > 0.9 there are two contour lines (for each flux value) approximately symmetric 
with respect to Z = 90° (see Figure la). On the other hand, for q<0.9 there is a single contour line (for 
each value of the flux) which varies much less with the longitude (see Figure lb). 

As we can see from the Table and the Figures, the predicted value for the 7-ray flux at high-galactic 



latitude is very close to that found by Dixon et al. (1998). This conclusion holds almost irrespectively of 



the flatness parameter. Moreover, the comparison of the overall shape of the contour lines in our Figures 



la and lb with the corresponding ones in Figure 3 of Dixon et al. ( 1998 ) entails that models with flatness 
parameter q < 0.8 are in better agreement with data, thereby implying that most likely the halo dark matter 
is not spherically distributed. 



10 As it would be exceedingly difficult to keep track of the dumpiness of the actual gas distribution, we assume that its 
density goes like the dark matter density - anyhow, the very low angular resolution of 7-ray detectors would not permit to 
distinguish between the two situations. 

As usual, we use the coordinate transformation x = — pcosfecosZ + i?o, V = — pcosfesinZ and z = psinb, where Rq = 8.5 
kpc is our galactocentric distance. 
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We remark that eq. @ yields $ 7 DM (> 0.1 GeV) ~ 5.9 x 10 6 7 s 1 cm 2 sr 1 at high-galactic 
latitude (for a spherical halo) . This value is roughly 40% of the diffuse extragalactic 7-ray emission of 



1.45 ± 0.05 x 10 5 7 s 1 cm 2 sr 1 found by the EGRET team (Sreekumar et al. 1998 ). So, our result 



supports the conclusion of Dixon et al. ( 1998] ) that the halo 7-ray emission is a relevant fraction of the 



isotropic diffuse flux also for E 1 > 0.1 GeV. 

Before closing this Letter, we would like to briefly address the crucial question whether the newly 
discovered 7-ray halo emission really calls for a dark matter source. For, one might suspect that a 
nonstandard inverse-Compton 7-ray production mechanism could explain the data (owing to the large 
uncertainties both in the electron hight scale and in the electron injection spectral index). However, this 
seems not to be the case. Basically, the inverse-Compton contour lines decrease much more rapidly than 
the observed ones. Hence, it would be impossible to explain in this manner the 7-ray flux found by Dixon 



et al. (1998) while still correctly accounting for the observed disk emission (Sreekumar et al. 1998|) . A more 



detailed account of this topic will be presented elsewhere. 

In conclusion, we feel that - in spite of the various uncertainties - the remarkably good agreement 
between theory and experiment makes our model for halo dark matter worth further consideration. In 
particular, the next generation of 7-ray satellites like AGILE and GLAST can test our prediction, thanks to 
the higher sensitivity and the better angular resolution. In this respect, it might be interesting to measure 
whether there is an enhancement in the 7-ray flux towards the nearby M31 galaxy, since we expect a similar 
mechanism for the 7-ray production to hold in its halo as well. 

The work of FDP is supported by an INFN grant. We would like to thank G. Bignami, P. Caraveo, D. 
Dixon, T. Gaisser, M. Gibilisco, G. Kanbach, T. Stanev, A. Strong and M. Tavani for useful discussions. 
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Table 1. Halo 7-ray intensity at high-galactic latitude for a spherical halo evaluated for R m in = 10 and 
15 kpc at energies above 0.1 GeV and 1 GeV for different values of the CR spectral index a (see eq. (ph). 



Rmin E 1 a $ 7 DM (fe = 90°) 



(kpc) 


(GeV) 




(7 cm 2 s 1 sr 1 ) 




10 


0.1 


2.45 
2.70 
3.00 


6.2 x 10 -6 
5.9 x 10~ 6 
4.9 x 10" 6 


10 


1 


2.45 
2.70 
3.00 


1.1 x 10 _B 
6.7 x 10~ 7 
3.3 x 10" 7 


15 


0.1 


2.45 
2.70 
3.00 


3.7 x 10~ 6 
3.5 x 10~ 6 
2.9 x 10" 6 


15 


1 


2.45 
2.70 
3.00 


6.5 x 10" 7 
4.0 x 10~ 7 
1.9 x 10" 7 



